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Abstract 

Recent data from the Tevatron has revealed that the produc- 
tion rate of prompt charmonium at large transverse momentum 
is orders of magnitude larger than the best theoretical predictions 
of a few years ago. These surprising results can be understood by 
taking into account two recent developments that have revolution- 
ized the theoretical description of heavy quarkonium production. 
The first is the realization that fragmentation must dominate at 
large transverse momentum, which implies that most charmo- 
nium in this kinematic region is produced by the hadronization 
of individual high-pT partons. The second is the development of a 
factorization formalism for quarkonium production based on non- 
relativistic QCD that allows the formation of charmonium from 
color-octet cc pairs to be treated systematically. This review sum- 
marizes these theoretical developments and their implications for 
quarkonium production in high energy colliders. 



1 INTRODUCTION 

The goal of high energy physics is to identify the elementary constituents 
of matter and to understand their fundamental interactions. Over the 
last twenty years, this endeavor has been extraordinarily successful. A 
gauge theory called the minimal Standard Model provides a satisfactory 
description of the strong, weak, and electromagnetic interactions of all 
the known elementary particles. There are very few discrepancies be- 
tween theory and experiment, and most of them are at the level of a few 
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standard deviations or less. However there is one process for which ex- 
perimental results have differed from theoretical predictions by orders of 
magnitude: the production of charmonium at large transverse momen- 
tum at the Tevatron. This dramatic conflict between experiment and 
theory presents a unique opportunity to make a significant step forward 
in our understanding of heavy quarkonium physics. 

Prior to 1993, the conventional wisdom on the production of char- 
monium in hadron collisions was based primarily on calculations in the 
color-singlet model. In this model, the production of a charmonium state 
is assumed to proceed through parton processes that produce a cc pair 
in a color-singlet state. It is the predictions of the color-singlet model 
at lowest order in as that disagree so dramatically with the Tevatron 
data. The data can, however, be explained by combining two recent 
theoretical developments in heavy quarkonium physics. The first is the 
realization that heavy quarkonium at large transverse momentum is 
produced primarily by fragmentation, the hadronization of individual 
high-pT partons. In the color-singlet model, fragmentation first con- 
tributes at higher order in as and thus was not taken into account in 
previous calculations. The second development is the realization that 
color- octet mechanisms, in which the cc pair is produced at short dis- 
tances in a color-octet state, sometimes dominate the production of 
charmonium. These mechanisms can be analyzed systematically using 
a new factorization formalism for heavy quarkonium production that is 
based on an effective field theory called nonrelativistic QCD (NRQCD). 
The Tevatron data can be explained by including a color-octet term in 
the fragmentation function for the formation of charmonium in a gluon 
jet. Evidence in support of this explanation has been accumulating, but 
a great deal of work, both experimental and theoretical, will be required 
before it can be regarded as conclusive. 

The purpose of this review is to summarize the recent theoretical 
developments in heavy quarkonium production. A brief review of these 
developments has been given previously by Mangano ^ . Our main focus 
will be on the production of quarkonium in high energy colliders, where 
both fragmentation and color-octet mechanisms are important. The 
color-singlet model is reviewed in Section 2. We summarize its predic- 
tions for charmonium production in pp colliders, which are in dramatic 
disagreement with Tevatron data. The production of quarkonium by 
fragmentation is reviewed in Section 3. We describe the factorization 
theorems of perturbative QCD that imply that fragmentation should 
dominate at large transverse momentum, and we illustrate them us- 
ing the electromagnetic production of the J/tp in decay. We then 
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discuss the predictions of the color-singlet model for fragmentation func- 
tions. Color-octet mechanisms for producing quarkonium are reviewed 
in Section 4. We first discuss the production of P-wave states, where 
a color-octet mechanism is required for perturbative consistency. We 
then describe the NRQCD factorization formalism, which implies that 
color-octet mechanisms must also contribute to the production of S-wave 
states. We also discuss its implications for fragmentation functions. In 
Section 5, we summarize applications of the recent theoretical devel- 
opments to the production of prompt charmonium, bottomonium, and 
the Be at the Tevatron and at LEP. We conclude in Section 6 with an 
outlook on the work that will be required to develop a comprehensive un- 
derstanding of heavy quarkonium production in high energy processes. 

2 THE PROBLEM OF CHARMONIUM 
PRODUCTION 

Before 1993, most predictions for charmonium production were based 
on the color-singlet model. In this Section, we review the color-singlet 
model and summarize its predictions at leading order in as for charmo- 
nium production at large transverse momentum in pp collisions. We then 
describe the experimental results from the Tevatron that have forced a 
reexamination of the problem of charmonium production. 

2.1 Color- singlet Model 

It is difficult to ascribe credit for the color-singlet model, since many 
physicists were involved in its early development. The decays of B 
mesons into charmonium states through the process 6 — > cc-\- s was first 
treated in the color-singlet model by DeGrand and Toussaint, by Wise, 
and by Kiihn, Nussinov, and Riickl ||^. The hadronic production of JM 
through the parton process gg ^ cc -\- g was calculated by Chang |g] . 
A thorough treatment of charmonium production in hadron collisions 
through 2^3 parton processes was later presented by Baier and Riickl 
Q. Guberina, Kiihn, Peccei, and Riickl Q applied the color-singlet 
model to charmonium production from the decay Z'^ — s- cc -I- 7. Berger 
and Jones Q calculated the rate for photoproduction of charmonium, 
and emphasized the kinematical restrictions on the applicability of the 
color-singlet model. This model was also applied to the inclusive pro- 
duction of charmonium in e'^e~ annihilation by Keung and by Kiihn 
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and Schneider |^] . A thorough review of the apphcations of the color- 
singlet model to heavy quarkonium production was recently given by 
Schuler @. 

An alternative model for quarkonium production called the color- 
evaporation model was developed around the same time In this 
model, it was assumed that all cc pairs with invariant mass between 2mc 
and the DD threshold 2m/) produce charmonium states. The fraction 
fn of the cc pairs that form a particular charmonium state H is assumed 
to be independent of the production process. This model is incapable of 
describing the variation of the production ratios for charmonium states 
between processes and as functions of kinematical variables. It will 
therefore not be considered further. 

To motivate the color-singict model, we can think of the production 
of charmonium as proceeding in two steps. The first step is the pro- 
duction of a cc pair, and the second step is the binding of the cc pair 
into a charmonium state. In order to have a significant probability of 
binding, the cc pair must be produced with relative momentum that, 
in the cc rest frame, is small compared to the mass mc of the charm 
quark. Otherwise, the c and c will fly apart and ultimately form D and 
D mesons. 

We first consider the production of the cc pair. Assuming that the c 
and c are not present in the initial state, any Feynman diagram for the 
production of a cc pair must involve virtual particles that are off their 
mass shells by amounts of order nic or larger. The part of the amplitude 
in which all internal lines are off-shell by amounts of order rUc or larger 
is called the short- distance part, and it is calculable using perturbation 
theory in as{mc). The parts of the amplitude in which the c and c lines 
are off-shell by amounts much less than mc can be considered part of 
the amplitude for the formation of the bound state. The short-distance 
part of the amplitude describes the production of a cc pair with a spatial 
separation that is of order l/rUc or smaller. This follows from the fact 
that the short-distance part is insensitive to changes in the relative 3- 
momentum of the c and c that are much less than mc- Since l/rUc 
is much smaller than the length scale associated with the charmonium 
wavefunction, the cc pair is essentially pointlike on that scale. Thus 
we need only consider the amplitude for a pointlike cc pair to bind to 
form a charmonium state. This amplitude will necessarily depend on 
the charmonium state H and on the quantum numbers of the cc pair. 

For any given charmonium state, the dominant Fock state consists of 
a color-singlet cc pair in a definite angular-momentum state. We denote 
the two possible color states of a cc pair by 1 for color-singlet and 8 for 
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color-octet. We use the spectroscopic notation ^^'^^Lj for the angular 
momentum state, where S, L, and J are the quantum number associ- 
ated with the total spin, the orbital angular momentum, and the total 
angular momentum, respectively. Thus the dominant Fock state for a 
charmonium state H is denoted |cc(]_, ^"^''"^Lj)) for some appropriate 
values of S, L, and J. For example, the dominant Fock state for the 
J/V; is \cc{l,^Si)), while for the Xcj, it is |cc(l, ^Pj)). 

The color-singlet model is a simple model for the amplitudes for a 
pointlike cc pair to form a charmonium state. If the dominant Fock 
state of the meson H is \cc{l, ^■^+^Lj)), the amplitude is assumed to be 
unless the pointlike cc pair is in a color-sing let ^■^+^Lj state. For this 
state, the amplitude can be expressed in terms of the L'th derivative of 
the radial wavefunction at the origin for the meson H. For example, the 
amplitudes for producing the states ij} and Xcj, plus some specific final 
state F , are assumed to have the forms 

A{^ + F) = A{cc{l,^Si)+F) R^{Q), (1) 

A{xcj + F) - A{cc(l,''Pj)+F)R'^^{Q). (2) 

The ^'s are amplitudes for producing color-singlet cc pairs with van- 
ishing relative momentum in the angular-momentum states indicated. 
The factor i?^(0) is the radial wavefunction at the origin for the ip, while 
R'^^ (0) is the derivative of the radial wavefunction at the origin for the 
Xc states. In the color-singlet model, it is assumed that the amplitudes 
A can be calculated using perturbative QCD and that all nonpcrtur- 
bative effects can be absorbed into the wavefunction factors. From (|l|) 
and (H), we deduce that the corresponding inclusive differential cross 
sections in the color-singlet model have the forms 

da{^ + X) = da(cc(l,35i)+X) |i?^(0)|2, (3) 
da{xcj + X) dd{cc{l,^Pj)+X)\R'^^{0)\^. (4) 

The color-singlet model has enormous predictive power. The cross 
section for producing a quarkonium state in any high energy process 
is predicted in terms of a single nonperturbative parameter for each 
orbital-angular-momentum multiplet. For example, the nonperturbative 
factor is i?^(0) for the S-wave states and rjc- It is -R^^(O) for the P- 
wave states XcOj Xci, Xc2, and he- Moreover, these parameters can be 
determined from decays of the charmonium states. For example, i?^(0) 
can be determined from the electronic width of the ^: 

r(^-.c+e-) « ^\Rm'- (5) 
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Thus the color-singlet model gives absolutely normalized predictions for 
the production rates of charmonium states in high energy collisions. 

In spite of its great predictive power, the color-singlet model is only 
a model. There are no theorems that guarantee that the amplitude 
factors in the simple way assumed in (]l|) and (H). In particular, it 
was never proven that higher order radiative corrections would respect 
the factored form. In addition, the model is clearly incomplete. For 
one thing, relativistic corrections, which take into account the relative 
velocity v of the quark and antiquark, are neglected. These corrections 
are probably not negligible for charmonium, since the average value of 

is only about 1/3. The color-singlet model also assumes that a cc 
pair produced in a color-octet state will never bind to form charmonium. 
This assumption must break down at some level, since a color-octet cc 
pair can make a nonperturbative transition to a color-singlet state by 
radiating a soft gluon. The most glaring evidence that the color-singlet 
model is incomplete comes from the presence of infrared divergences 
in the cross sections for P-wave states. This problem and its solution 
will be discussed in Section 4. For the moment, we simply note that 
the infrared divergence violates the factorization assumption implicit in 
(^). It implies that the cross section da is sensitive to small momentum 
scales, so it cannot be calculated reliably using perturbative QCD. 

2.2 Expectations for Charmonium at Large pt 

Most calculations of charmonium production prior to 1993 were based 
on two crucial assumptions. The first was that the amplitude for the 
formation of a charmonium state from a cc pair was accurately de- 
scribed by the color-singlet model. The second was that the dominant 
production processes for color-singlet cc pairs were the Feynman dia- 
grams that were lowest order in a^. The first thorough treatment of 
the problem of charmonium production at large transverse momentum 
in hadron-hadron collisions was given by Baier and Riickl in 1983 
In subsequent calculations [0, the contributions from i?-meson decay 
into charmonium states were included. The results of these calculations 
are summarized below. 

We first introduce some terminology for describing charmonium pro- 
duction in high energy colliders. Charmonium is called prompt if the 
point at which the charmonium state is produced and the collision 
point of the colliding beams can not be resolved using a vertex detec- 
tor. Prompt charmonium is produced by QCD production mechanisms. 
The decay of 6-hadrons produces charmonium that is not prompt. If 



8 



BRAATEN, FLEMING & YUAN 



a 6-hadron is produced with large transverse momentum, it will travel 
a significant distance before decaying weakly. For a 6-hadron with pt 
around 10 GcV, the displacement between the collision point and the 
secondary vertex where the charmonium is produced is typically a frac- 
tion of a millimeter. 

A charmonium state that is prompt, but does not come from the 
decay of a higher charmonium state, is called direct. For example, the 
ij) is produced in the decays of x^o, Xci) Xc2, and V' with branching 
fractions of approximately 0.7%, 27%, 14%, and 57%, respectively. Thus 
the prompt ■0 signal includes direct V''s and contributions from direct 
Xcj and direct tf}' . 

A charmonium state with large momentum is called isolated if there 
are no other hadrons whose momentum is nearly collinear. Charmonium 
that is produced from the decay of a 6-hadron with large transverse 
momentum is never isolated, because the remnant hadrons from the 
decay of the 6-hadron will have momentum that is nearly collinear. 

We now consider the mechanisms for the production of charmonium 
with large px in pp collisions. First we focus on charmonimn production 
from the weak decay of 6-hadrons. Note that the decay of 6-hadrons 
with large pr produces charmonium that is neither prompt nor isolated. 
The dominant mechanism for producing 6-hadrons at large p^ is the 
gluon fusion process gg bb, followed by the hadronization of the 
b or b. The inclusive branching fractions for any 6-hadron to decay 
into charmonium states are presumably close to those for B mesons. 
The inclusive branching fractions for S-meson decays are approximately 
1.3% for V, 1% for Xci, and 0.5% for V'- 

The color-singlct model gives predictions for the production rate 
of charmonium from QCD mechanisms. If a cc pair is produced with 
large Pt through parton collisions, there must be a recoiling parton to 
balance the transverse momentum. Thus the production mechanisms 
that are of leading order in as are 2^3 parton processes. Note that 
the charmonium states produced by these processes are prompt and 
isolated. In the color-singlet model, the only 2^3 process that can 
produce ip or ip' is gg ^ cc + g. The 2^3 processes that produce 
Xcj's with large pt are gg ^ cc + g, gq ^ cc + q, gq ^ cc + q, and 
qq cc + g. The cross sections for these 2^3 processes arc all of order 
al- At large px, the parton differential cross sections da/dp^ scale like 
l/p|, for the XcJ states and like l/p% for ij) and if)'. Thus the parton 
cross sections for ^ and ip' are suppressed relative to those for XcJ by a 
factor of mg/p|n at large pr- 

We now summarize the predictions of the color-singlet model at lead- 
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ing order in for inclusive production in pp collisions |1C| . The 
contribution from 6-hadron decay falls off the most slowly with px and 
it was predicted to dominate at the Tevatron for px > 7 GeV. Of the 
prompt production mechanisms for ■(/;, the most important was found to 
be the decay of direct Xci ■ The contributions from direct Xc2 and direct 
tp were down by factors of about 4 and 18 at px — 10 GeV. Thus the 
conventional wisdom before 1993 was that iJj production at large pr at 
the Tevatron should be dominated by 6-hadron decay, with decays of 
direct Xci and direct Xc2 being the only other important mechanisms. 

In the case of tp' , the decay of 6-hadrons was predicted to be the only 
important production mechanism at large px- It should be emphasized 
that these predictions were based on the color-singlet model and on the 
additional assumption that the dominant parton processes were of order 

2.3 Prompt Charmonium at the Tevatron 

The first substantial data on the production of charmonium in pp col- 
lisions came from the SppS operating at a center-of-niass energy of 630 
GeV. Data collected by the UAl collaboration contained indications 
of deviations from the predictions of the color-singlet model jl^ . The 
Tevatron, which operates at the significantly higher energy of 1.8 TeV, 
has provided an opportunity to investigate charmonium production in 
much greater detail. The CDF collaboration has accumulated large sam- 
ples of data on the production of '0, XcJi and ip' at the Tevatron. Though 
analysis of the data from the 1988-1989 collider run was hampered by 
the inability to separate prompt charmonium from charmonium that is 
produced by the decays of 6-hadrons, discrepancies between the data 
and the predictions of the color-singlet model were evident ||l^. How- 
ever, assuming the conventional wisdom that ip production at large pr 
should be dominated by the decays of 5-hadrons and direct Xc's, CDF 
extracted a value for the b quark cross section from their data on and 
Xc production Q. This value was shown to be too large by about a 
factor of 2 in the subsequent collider run iQ . 

Before the 1992-93 run of the Tevatron, CDF installed a silicon ver- 
tex detector that can resolve secondary vertices separated by distances 
greater than about 10 /xm from the pp collision point ]l5| . Prompt char- 
monium is produced essentially at the collision point, while charmonium 
from the decay of 6-hadrons is produced at a secondary vertex which 
is typically hundreds of microns away. Thus the silicon vertex detector 
can be used to separate prompt charmonium from charmonium that is 
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produced in 5-liadron decays. 

In the charmonium data sample collected during the 1992-93 run, 
it was found that only 20% of the ip's and only 23% of the tp'^s come 
from the decay of 6-hadrons ||l^. Furthermore, only 32% of the ^''s 
come from Xc decay [0. Thus, the majority of the ip^s and 'ip"s must 
be produced by other mechanisms, in dramatic contradiction to the 
conventional wisdom. Furthermore, the fraction of ip and tp' from b- 
hadron decay does not increase appreciably with increasing pj^ , contrary 
to the predictions of the color-singlet model. 

An excellent review of the results on quarkonium production from 
the Tevatron, and also from other experiments, has been presented by 
Sansoni |l^. CDF has measured the cross sections for prompt ?/' and 
prompt ip' production as a function of px | |T^ . The prompt ip signal has 
been resolved into those ip's that come from Xc decays and those that 
do not ]l7t . In addition, the ratio of the production rates of Xci and 
Xc2 has been measured. The most recent CDF data on the production 
of prompt i/i's that do not come from Xc decay, prompt ijj^s from Xc 
decay, and prompt V'''^, are shown in Figures 1, 2, and 3, respectively. 
The prompt ip^s that do not come from Xc decay include direct V's 
and ip^s from the decay of direct ip'^s. The predictions of the color- 
singlet model at leading order in as are shown as dashed lines. These 
predictions fall several orders of magnitude below the data at large pT- 
Thus the predictions of the color-singiet model at lowest order in as fail 
dramatically when confronted with the data on charmonium production 
from the Tevatron. 

3 FRAGMENTATION 

The first major conceptual advance in the recent revolution in heavy 
quarkonium production was the realization that fragmentation dom- 
inates at sufficiently large transverse momentum. Fragmentation is 
the formation of a hadron within a jet produced by a parton (quark, 
anti-quark, or gluon) with large transverse momentum. As the parton 
emerges from the collision point, it radiates gluons and other partons, 
most of which are almost coUinear. The partons ultimately coalesce into 
the hadrons that make up the jet. In the case of production of charmo- 
nium by fragmentation, the "jet" containing the charmonium state may 
not qualify as a jet by conventional experimental definitions. For exam- 
ple, if most of the momentum of the "jet" is carried by the charmonium 
state, it may not satisfy a jet criterion that requires a specified number 
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of tracks above a certain momentum threshold. 

The word fragmentation is sometimes used to describe the coales- 
cence of partons into hadrons, whether or not these partons make up 
a jet. We prefer to use the word hadronization to describe this general 
process. The word fragmentation will be reserved specifically for the 
formation of a hadron within the jet produced by a high-pT parton. 
Thus fragmentation involves the hadronization of the partons in the jet, 
but hadronization also occurs in low energy processes that have nothing 
to do with jets. Fragmentation is a useful concept because the proba- 
bility for the formation of a hadron within a jet is universal, i.e. it is 
independent of the process that produces the parton that initiates the 



In this Section, we introduce the factorization theorems of pertur- 
bative QCD that guarantee that inclusive hadron production at large 
transverse momentum is dominated by fragmentation. We illustrate the 
factorization theorems using the simple example of the electromagnetic 
production of tp in decay. Finally we discuss the color-singlet model 
predictions for the fragmentation functions of heavy quarkonium. 

3.1 Factorization Theorems of Perturbative QCD 

One of the classic factorization theorems of perturbative QCD guar- 
antees that inclusive hadron production in e+e~ annihilation at suffi- 
ciently large energies is dominated by fragmentation. After describing 
this factorization theorem, we discuss its extension to inclusive hadron 
production at large transverse momentum in hadron-hadron collisions. 

We consider the inclusive production of a hadron H with energy E in 
e^e^ annihilation at large center-of-mass energy ^/s. We are interested 
in the cross section in the scaling limit in which E, ^/s — > oo with E/^/s 
held fixed. The production of the hadron also involves lower momentum 
scales, such as the hadron mass and the scale Aqcd associated with 
nonperturbative effects in QCD. The factorization theorem states that 
the cross section in the scaling limit has the form [p^ 



jet. 



d(T{e^e 



H{E)+X) 




(6) 



where the sum is over parton types i and the integral is over the longi- 
tudinal momentum fraction z of the hadron H relative to the parton i. 
In (ph, d<7 is the differential cross section for producing a parton i with 
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total energy E/ z. This cross section is only sensitive to momenta on the 
order of E, so it can be calculated using QCD perturbation theory. The 
effects of lower momentum scales can be systematically factored into 
the functions ZJ^^/f (z, /i), which are called fragmentation functions or 
parton decay functions. The fragmentation function Di^niz, i^) gives 
the probability that the jet initiated by parton i will include a hadron 
H carrying a fraction z of the jet momentum. The factorization the- 
orem holds to all orders in perturbation theory. For a light hadron H 
whose mass is on the order of Aqcd, all corrections to the factorization 
formula (||) fall like powers of Aqcd /E. 

The essential ingredient in the proof of the factorization theorem 
was given by Collins and Sterman . They demonstrated that a dia- 
gram that contributes to the inclusive cross section in the scaling limit 
can be separated into a hard-scattering subdiagram that produces hard 
partons, jet-like subdiagrams for each of the hard partons, and a soft 
part. The soft part includes soft gluon lines that can couple to any of 
the jet-like subdiagrams. After summing over all possible connections of 
the soft gluons, one finds that the effects of the soft parts cancel, leaving 
a factored form for this contribution to the inclusive cross section. In 
the case of inclusive hadron production, these ideas were used by Curci, 
Furmanski, and Petronzio and by Collins and Soper Q to provide field 
theoretic definitions of the fragmentation functions. 

The factorization theorem (|6|) requires the introduction of an arbi- 
trary scale ^ that separates the large momentum scale E from the lower 
momentum scales. The parton cross sections and the fragmentation 
functions depend on the arbitrary scale n in such a way that the cross 
section is independent of /i. The /x-dependence of the fragmentation 
functions is given by an evolution equation of the form 



The kernel Pi^j {x, ^) describes the splitting of a parton i into a parton j 
with momentum fraction x, and is calculable as a perturbation series in 
as(/x). At leading order in a^, these kernels are identical to the Altarelli- 
Parisi functions that govern the evolution of parton distributions. 

The factorization theorem for inclusive hadron production in e'^e~ 
annihilation can be generalized to other high energy processes. The real 
power of these factorization theorems lies in the fact that the fragmenta- 
tion functions are universal, i.e. they are independent of the process that 
produces the fragmenting partons. Thus, if the fragmentation functions 
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for a hadron H are determined from e^e~ annihilation data, they can 
be used to predict the production rate of the hadron H in jets produced 
by other high energy processes. 

One of the important generahzations of the factorization theorem 
is to inclusive hadron production at large transverse momentum pT in 
hadron- hadron collisions. The proof of the factorization theorem for 
this process is more difficult, because there are two hadrons in the ini- 
tial state. In fact, such a proof has actually been carried out only for 
the simpler case of the Drell-Yan process for creating a muon pair p^ . 
However there are no apparent obstacles to extending this proof to the 
case of inclusive hadron production at large px jl^ . The resulting fac- 
torization formula is 

da{AB ^ H{pt)+X) = f dxi fj/A{xi) f dx2 fk/sixi) 

X I dz daijk ^ i{pT/z) + X) A^ff(^)- (8) 
Jo 

Long-distance effects can be factored into the fragmentation functions 
and into the parton distributions fj/Ai^i) and fk/B{x2) for hadrons A 
and B. The restriction to large px is necessary in order that the jets 
consisting of the remnants of hadrons A and B after the hard scatter- 
ing have large momentum relative to the jet containing hadron H . The 
formula (||) should hold to all orders in perturbation theory, with correc- 
tions falling like powers of ^qcd/pt- Implicit in (H) are three arbitrary 
scales: the factorization scale ^p, which cancels between the parton dis- 
tributions and c?(T, the fragmentation scale fifrag, which cancels between 
da and the fragmentation functions, and the renormalization scale fj-n 
for the running coupling constant, which appears in da. In low-order 
calculations, Hhag , A*f, A*;? should all be chosen on the order ofpr/z, 
the transverse momentum of the fragmenting parton. 

In 1993, Braaten and Yuan pointed out that the factorization 
theorems for inclusive hadron production must apply to heavy quarko- 
nium as well as to light hadrons. The only difference is that the lead- 
ing corrections fall as powers of uiq/pt, where mg is the heavy quark 
mass, instead of powers of J^qcd/pt- Therefore the dominant produc- 
tion mechanism for charmonium at pt 3> Wc must be fragmentation. In 
retrospect, this statement may seem obvious, but fragmentation contri- 
butions were not included in any of the previous calculations of charmo- 
nium production in pp collisions summarized in Section 2.2. This can 
easily be seen from the pt dependence of the parton cross sections. In 
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the factorization formula the only momentum scale that the par- 
ton cross section da can depend on is pt- Therefore, by dimensional 
analysis, da/dp"^ scales like at large px- The color-singlet model 

cross sections at leading order in fall off much more rapidly with px- 
da/dp^ scales like 1/p^ for -0 and i>' and like l/pf^ for Xcj- 

3.2 Electromagnetic Production of i/j in Decay 

To illustrate the factorization theorem, we discuss the production of ip in 
decay through electromagnetic interactions. This example illustrates 
two important points. First, the production process that is lowest order 
in the coupling constant does not necessarily dominate in the asymptotic 
region. Second, the terms that do dominate in the asymptotic region can 
be factored into cross sections for producing partons and fragmentation 
functions. 

The electromagnetic production of ijj is particularly simple, because 
the QCD interaction can be treated nonperturbatively by expressing the 
production amplitude in terms of a matrix element of the electromag- 
netic current between the QCD vacuum and a Tp: 

(V|fc7^c|0) = g^My, (9) 

where is the polarization vector for the ip and is dimensionless. 
The value of 5^ can be determined from the decay rate for ^l) e~^e~: 

r(V^e+e-) = —aHl^gl. (10) 

This gives the value g^ = 0.008. 

The various contributions to electromagnetic ip production in de- 
cay can be calculated using perturbation theory in the electromagnetic 
coupling constant a. Since this coupling constant is small, we might 
naively expect the dominant production process to be the one that is 
lowest order in a. We must remember, however, that there is also an- 
other small dimensionless parameter in this problem, namely M^/Mz- 
The relative importance of a production process will be determined both 
by its order in a and by how it scales with /Mz- We will see that the 
leading-order process is suppressed by AI^/M'^, and that the dominant 
production process is actually higher order in a. 

The production process for ^p that is lowest order in a is Z'^ ^ ip^, 
which proceeds at order through the process Z^ ^ cc + 7. For the 
cc to form a ip), the relative momentum of the c and c in the ip rest 
frame must be small compared to nic- This is in turn small compared 
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to the momentum of the ip, which is approximately AIz /2. If the relative 
momentum of the c and c is neglected, the amplitude is proportional to 
the matrix element {ip\cj^c\Q), and it therefore can be expressed in terms 
of g^. The branching fraction for Z° — > ^pJ is 5.2 x 10-8 §. In the limit 
<C Mz, the branching fraction is proportional to ag'^M^/M^, so 
it vanishes in the scaling limit. The suppression factor M^/M'^ reflects 
the fact that the c or c must receive a momentum kick of order Mz from 
the photon in order for the c and c to be produced with small enough 
relative momentum to form a bound state. 

At next-to-leading order in a, the ip can be produced electromagnet- 
ically via the decay Z° iIj + £'^£~, which proceeds at order a"^ through 
the process cc + t^i^ . The branching fraction for Z° ip + £~^£~ 

is 7.5 X 10"^ |2^. Thus the decay rate for this order-a'^ process is 
an order of magnitude larger than the order-a^ process Z'^ — s- tpj. The 
reason it is larger is that the factor of M^/M'^ in the decay rate for 
Z'^ ipj provides larger suppression than the extra factor of a in the 
decay rate for — -0 -I- £'^£~ . 

We introduce a scaling limit for the production of -ip with energy E 
in Z° decay. This hmit is E, Mz oo with E/Mz and held fixed. 
In the scaling limit, the differential decay rate for the electromagnetic 
production of tp satisfies a factorization theorem that is analogous to 
that in (|): 

dr(z° -> i}{E) + x) 

= y2 f dzdfiZ'' ^i{E/z)+X,fi) D,^^{z,fi), (11) 

^ Jo 

where the sum over partons i includes photons and positive and negative 
leptons. The factor dT is the inclusive rate for decay into a parton i 
with energy E/z. The factor Di^^{z) is a fragmentation function that 
gives the probability for an electromagnetic jet initiated by the parton 
i to include a ■0 that carries a fraction z of the jet momentum. 

Fleming Eq] pointed out that, since the differential decay rate for 
the process Z^^ ip + £'^£~ does not vanish in the scaling limit, it must 
be expressible in the form (^ij). He found that, at leading order in a 
and in the scaling limit, it can be written as 

^{Z'' ^^{z) +£+£-) = 2r{Z°-^£+£-)De^^{z,fi) 
az 

+ ^{Z''^^{z)+£+£-,^^)P^^^, (12) 
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where z = 2E/Mz- The photon fragmentation probabiHty is P- 



Anag'^ = 7 x 10 The lepton fragmentation function is 



Di>^^{z,fi) = 2a 



iz-ir + 



1 Zf,' 



(13) 



where /i is an arbitrary factorization scale that separates the scales 
Mz and M^. The /i-dependence of the fragmentation function can- 
cels against that of dT/dz in (|l^). Note that all the effects of QCD, 
both perturbative and nonperturbative, are absorbed into the factor g'^ 
in Py^^ and in (p^). In the decay rate for ip + i^l^ , the absence 

of kinematic suppression factors like M^/M"^ can be attributed to the 
fact that the cc pair that forms the ip is produced at a momentum scale 
of order rric rather than Mz- The large momentum scale Mz only enters 
in the decay of the Zo into the partons £^ , and 7. 



3.3 Fragmentation Functions in the Color-singlet Model 

In 1993, Braaten and Yuan pointed out that the color-singlet model 
gives predictions for the fragmentation functions for the formation of 
heavy quarkonium in quark and gluon jets |2^ . They calculated the 
fragmentation functions explicitly for 5 — > t^c, V' to leading order in a^. 
For example, the fragmentation function for g rjc is calculated from 
the parton process g cc + g and is of order aj. : 



Dg^^^{z,2mc) 

^ al{2m,)\R^M? 
247rm3 



[3z-2z2 + 2(l-z)log(l-z)] . (14) 



This is the fragmentation function at an initial scale of order 2mc. It 
can be evolved up to a higher scale /i by using the evolution equations 
(0) . The fragmentation function for g — > i/) is calculated from the parton 
process g — > cc -I- gg and is of order a?. The fragmentation functions 
for c ric,''P [0 ^^'^ ^ ^ Bc,B* [|8| were subsequently calculated 
by Braaten, Cheung, and Yuan. The fragmentation functions for b — > 
Bc,B* had actually been calculated earlier by Chang and Chen p9|] , 
although they did not note the universality of these functions. 

The fragmentation functions for P-wave quarkonium states have also 
been calculated in the color-singlet model to leading order in a^. These 



fragmentation functions include g — > XcJ 31 1, c — > Xcj, and b —^ 
P-wave Be |3^, Some of the fragmentation functions for D- 

wave states have also been calculated p^, |3^ . The color-singlet model 
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fragmentation functions for g — > Xcj are calculated from the parton 
process g —f cc + g and are logarithmically infrared-divergent at leading 
order in a^. The solution to this problem will be discussed in Section 
4. It should be noted that there is a discrepancy between the two 
calculations of the fragmentation functions for g XcJ that has not yet 
been resolved 

The color-singlet model fragmentation functions can be resolved into 
the contributions from each of the possible values of the helicity h of the 
charmonium state. The fragmentation functions for c — > V't, where 
ipL and ipT represent the longitudinal {h — 0) and transverse = 1) 
polarization components of the ip, were first calculated by Chen |32] and 
by Falk et al Jsi ] . The fragmentation functions for g XcJ [ 51 , |3^ , 
b ^ B* |39| , and b — > P-wave be states |3^ have also been 
resolved into the contributions from the individual helicitics. 

In the method developed by Braaten and Yuan, the fragmentation 
functions were extracted by taking the scaling limit of color-singlet 
model cross sections and expressing them in the factored form demanded 
by the QCD factorization theorems. This method would be rather cum- 
bersome beyond leading order in as. The fragmentation functions can 
also be calculated directly from the field theoretic definitions . 
This method provides a significant advantage for calculating fragmen- 
tation functions beyond leading order in a^. Higher order corrections 
may be particularly important if the leading order calculation gives a 
soft fragmentation function that is small for z near 1 H]. 



4 COLOR-OCTET MECHANISMS 

The second major conceptual advance in the recent revolution in heavy 
quarkonium production is the realization that color-octet mechanisms 
can be important. Contrary to the basic assumption of the color-singlet 
model, a cc pair that is produced in a color-octet state can bind to form 
charmonium. While a cc pair that forms charmonium must ultimately 
be in a color-singlet state, the short- distance part of the process can 
involve the production of a color-octet cc pair. 

In this Section, we first discuss the case of P-wave states, where 
perturbative consistency requires that a color-octet term be added to the 
color-singlet model cross section. We then discuss a general factorization 
formalism for quarkonium that is based on NRQCD. This formalism 
can be used to factor production cross sections into short-distance parts 
that can be computed using perturbative QCD and nonperturbative 
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NRQCD matrix elements. Color-octet mechanisms for S-wave states 
arise naturally in this framework. Finally, we discuss the implications 
of this formalism for fragmentation functions. 

4.1 Color-octet Mechanism for P-waves 

The most glaring evidence that the color-singlet model is incomplete 
comes from the presence of infrared divergences in the production cross 
sections for P-wave states. There are analogous divergences in the anni- 
hilation decay rates for P-wave states, which were discovered by Barbieri 
and collaborators as early as 1976 j42). They found that the decay rate 
for XcJ ~* 995 depends logarithmically on the minimum energy of the 
final state gluon. In phenomenological treatments of the annihilation 
decay rates of P-wave states, the infrared singularity was avoided by 
imposing an ad hoc infrared cutoff on the gluon energy. This cutoff was 
sometimes identified with the binding energy of the charmonium state, 
but without any good justification. 

The presence of infrared divergences in the color-singlet model is 
not a problem if the divergences can be absorbed into the nonperturba- 
tive factor in the cross section. For example, there are linear infrared 
divergences associated with the Coulomb singularity in radiative cor- 
rections to S-wave cross sections, but they can be factored into |i?(0)p. 
Similarly, the radiative corrections for P-waves include a linear infrared 
divergence that can be factored into |i?'(0)p. However, the structure 
of the logarithmic infrared divergences for P-waves is such that they 
cannot be factored into |i?'(0)p. 

The solution to this problem was given by Bodwin, Braaten, and 
Lepage in 1992 |Q. They noted that the infrared divergence in the 
annihilation rate for Xcj 99.9 arises when the final state gluon becomes 
soft. After radiating the soft gluon from either the c or the c in the 
color-singlet ^Pj bound state, the cc pair is in a color-octet ^Si state. 
It then annihilates through the process cc — > qq. In this region of phase 
space, the short-distance part of the decay amplitude is therefore the 
annihilation of a cc pair in a color-octet ^Si state. Thus, in addition 
to the conventional term in the color-singlet model, the factorization 
formula for the decay rate of the XcJ into light hadrons must include a 
second term. In the color-singlet model term, the short-distance factor 
is the annihilation rate of a cc pair in a color-singlet ^Pj state and the 
long-distance nonperturbative factor is proportional to \R'^^ (0)P- In the 
second term, the short-distance factor is the annihilation rate of a cc pair 
in a color-octet ^Si_ state and the long-distance factor is the probability 
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for the Xcj to contain a pointlike cc pair in a color-octet Si state. The 
color-octet term can be interpreted as a contribution to the decay rate 
from the ccg component of the Xcj wavefunction. 

The solution to the problem of infrared divergences in the production 
cross sections for the XcJ states is similar p4[ . The infrared divergence 
arises from the radiation of a soft gluon from the c or c that form the 
color-singlet ^Pj bound state. Before the radiation, the cc pair is in 
a color-octet '^S'l state. The short-distance part of the process in this 
region of phase space is the production of a cc pair with small relative 
momentum in a color-octet ^Si state. The factorization formula for 
P-wave cross sections must therefore include two terms: 

da{xcJ + X) = dd{c-c{h^Pj)+X)\R'^^m^ 

+ (2J+1) dd{cc{%,^Si)+X) {Of). (15) 

The first term is the conventional term of the color-singlet model. In 
this term, the short-distance factor is the cross section for producing a 
cc pair with small relative momentum in a color-singlet ^Pj state. The 
second term in ( p^ ) represents a contribution to the cross section from 
a color- octet mechanism. The short-distance factor in this term is the 
cross section for producing a cc pair with small relative momentum in 
a color-octet ^Si state. The long-distance nonperturbative factor (Of) 
is proportional to the probability for a pointlike cc pair in a color-octet 
■^5*1 state to bind and form a Xc- The factorization formula (|l5| ) has 
been applied to Xd production in S-meson decays in T decays 
^5|] , and in photoproduction p6|] . 

Both the color-singlet factorization formula (||) for S-waves and the 
factorization formula (^|) for P-waves hold to all orders in perturbation 
theory in the nonrelativistic limit. This is the limit in which the typical 
velocity v of the heavy quark in the bound state goes to 0. Unlike 
the case of S-waves, a color-octet production mechanism is required by 
perturbative consistency in the case of P-wave states. 

This can be understood from the NRQCD factorization approach 
discussed in Section 4.2. This approach indicates that the color-singlet 
term for P-waves is suppressed by a factor of compared to S-waves. 
It is this suppression of the color-singlet term that makes it necessary 
to include the color-octet term in the case of P-waves. 

4.2 NRQCD Factorization 

The simple factorization formulas (^) for S-waves and (|l|) for P-waves 
are correct in the nonrelativistic limit, i.e. in the limit w — > 0, where v 
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is the typical relative velocity of the c and c in charmonium. Potential 
model calculations indicate that is only about i for charmonium and 
about for bottomonium. Thus the neglect of relativistic corrections 
suppressed by powers of introduces a systematic theoretical error 
that may not be negligible. For specific processes, the problem may be 
worse than suggested by the magnitude of . If there are other small 
parameters in the problem, such as as(m^ or kinematical parameters 
such as m^/p^, the terms that are of leading order in may be sup- 
pressed by these other parameters. If this suppression is large enough, 
terms that are subleading in may actually dominate. 

The key to understanding the structure of the relativistic correc- 
tions to quarkonium production is to unravel the momentum scales in 
the problem. If there is a very large momentum scale that is set by the 
kinematics of the production process, such as Mz in the case of de- 
cay or pt in the case of production at large transverse momentum in pp 
collisions, this scale can be removed from the problem using the factor- 
ization theorems of perturbative QCD discussed in Section 3. The next 
largest momentum scale is the heavy quark mass ttt-c, which is certainly 
important in the production of a cc pair. The actual formation of the 
bound state from the cc pair involves smaller momentum scales, includ- 
ing the typical momentum rricV of the heavy quark in quarkonium, its 
typical kinetic energy rricv'^ , and the scale Aqcd of generic nonpertur- 
bative effects in QCD. In charmonium and bottomonium, the mass mg 
is large enough that effects associated with this scale should be calcula- 
ble using perturbation theory in the running coupling constant as{mQ). 
Nonperturbative effects become important at smaller momentum scales. 

A powerful tool for isolating the effects of the scale uiq is nonrela- 
tivistic QCD (NRQCD), an effective field theory developed by Caswell 
and Lepage |47j. NRQCD is a formulation of QCD in which the heavy 
quark and antiquark are treated nonrelativistically. They are described 
by a Schroedinger field theory with separate 2-component spinor fields, 
rather than a single 4-component Dirac field. The gluons and light 
quarks are described by the relativistic lagrangian for ordinary QCD. 

NRQCD can accurately describe the physics of heavy quarkonium at 
energies that are much less than mg above the rest mass. This follows 
from the fact that virtual states with energies of order mg or larger have 
lifetimes of order 1/mg or smaller. In this time, quanta can propagate 
over distances that are at most of order l/mg. Since this distance is 
much smaller than the size of the quarkonium bound state, the effects of 
these virtual states can be taken into account through local interaction 
terms in the NRQCD lagrangian. 
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The inverse of the mean radius of 

a quarkonium state is a dynamicaUy generated momentum scale that 
is much smaller than its rest mass. The ratio of these momentum scales 
defines a small parameter v which can be identified with the typical ve- 
locity of the heavy quark in the quarkonium state. This small parameter 
V, which is defined nonperturbatively, can be exploited to organize cal- 
culations of heavy quarkonium observables into expansions in powers of 
v^. NRQCD is equivalent to full QCD in the sense that the parameters 
in the NRQCD lagrangian can be tuned so that its predictions agree 
with those of full QCD to any desired order in |Q. 

Bodwin, Braaten, and Lepage have recently developed a theo- 
retical framework for inclusive quarkonium production that allows rel- 
ativistic corrections to be included to any desired order in w^. This 
framework provides a factorization formula for inclusive cross sections 
in which all effects of the scale TOc are separated from the effects of lower 
momentum scales, including rricV. The derivation of this factorization 
formula involves two steps. This first step is the topological factorization 
of diagrams that contribute to the cross section. This step is similar to 
the derivations of the factorization theorems of perturbative QCD that 
were discussed in section 3.1. For simplicity, we consider the case of 
charmonium production in e'^e~ annihilation at a center-of-mass en- 
ergy that is significantly larger than 2mc. We define a scaling limit 
by -^i, rric — > oo with y/s /rric fixed. A diagram that contributes to the 
cross section in this limit can be separated into a hard-scattering sub- 
diagram that produces hard partons and a cc pair with small relative 
momentum, a jet-like subdiagram for each of the hard partons, a subdi- 
agram that involves the c and c, and a soft part. The soft part includes 
soft gluons that can couple to the jet-like subdiagrams and to the cc 
subdiagram. After summing over all possible connections of the soft 
gluons, the effects of the soft gluons that are connected to the jet-like 
subdiagrams cancel, leaving a factored form for this contribution to the 
cross section. The derivation of topological factorization requires the 
hard partons to all have large momentum relative to the cc pair. This 
can presumably be generalized to processes in which there are hadrons 
in the initial state, provided that the charmonium is produced with 
transverse momentum large compared to Age d ■ 

After the topological factorization of the diagram, the hard-scatter- 
ing amplitude contains all effects of the scale rric, but it also depends on 
the relative momentum p of the cc pair. To complete the derivation of 
the factorization formula, the dependence on p must be removed from 
the hard-scattering amplitude, so that all effects of the scale rricV reside 
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in the cc subdiagram. This could be accomphshed simply by Taylor- 
expanding the hard-scattering amplitude in powers of p, if it were not for 
the fact that this generates ultraviolet divergences in the cc subdiagram. 
This is where the power of NRQCD becomes evident. This effective field 
theory can be used to systematically unravel the scales rric and mcU, so 
that all effects of the scale rricV are contained in the cc subdiagram. 
Essentially, the cc subdiagram is defined to include all the parts of the 
diagram that are reproduced by NRQCD. Those parts involving the cc 
pair that are not reproduced by NRQCD must involve the scale nic, and 
therefore can be included in the hard-scattering subdiagram. The result 
is a factorization formula for the inclusive cross section for producing a 
quarkonium state H that holds to all orders in a^. It has the form 

dcr{H + X) = ^dd{cc{n) + X) {O^), (16) 

n 

where da is the inclusive cross section for producing a cc pair in a color 
and angular-momentum state labeled by n and having vanishing rela- 
tive momentum. The parton cross sections da involve only momenta of 
order ttLc or larger, and therefore they can be calculated as perturbation 
expansions in as{mc). Since da is insensitive to relative momenta that 
are much smaller than mc, it describes the production of a cc pair with 
separation less than or of order 1/toc. This separation is essentially 
pointlike on the scale of a charmonium wavefunction, which is of order 
l/{mcv). Thus the nonperturbative long-distance factor {O^) is pro- 
portional to the probability for a pointlike cc pair in the state n to form 
the bound state H. Note that the only dependence on the quarkonium 
state H in (^6|) resides in the factor (O^). 

The factorization formula ( |l6|) holds only for inclusive cross sections. 
The matrix elements (O^) are proportional to the probabilities for the 
formation of the charmonium state H, plus light hadrons whose total 
energy in the H rest frame is of order rricV^. Summation over these 
light hadronic states is essential for the separation of short-distance and 
long-distance effects. The methods required to derive the factorization 
formula break down for exclusive cross sections. The derivation also 
implies that factorization does not hold at the amplitude level, contrary 
to the basic assumption of the color-singlet model. 

The factors {O^) in (|l6|) can be expressed as vacuum matrix ele- 
ments of 4-quark operators in NRQCD [Q. The operator creates a 
pointlike cc pair in the state n, projects onto states that in the asymp- 
totic future include the quarkonium state H , and finally annihilates the 
cc pair at the creation point, again in the state n. The matrix elements 
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that play the most important roles in quarkonium production were de- 
noted (Of and (Of in Ref. Q. The operator 
Oni^^'^^Lj) creates and annihilates a cc pair in the state cc(n, 
If the dominant Fock state of the meson H is \cc{l, "^^^^Lj)), the matrix 
element {O^ {^^^^ Lj)) can be related to the radial wavefunction. For 
example, up to corrections of relative order u"*, we have 

{Oti'S,)) = f |i?^(0)p, (17) 
zvr 

{Or'^Pj)) = (2J + l)|-|i?^jO)r (18) 

If we keep only these terms in the factorization formula (p^), we re- 
cover the formulas (^ and (||) of the color-singlet model. NRQCD has 
approximate heavy-quark spin symmetry, and this implies relations be- 
tween color-octet matrix elements. For example, up to corrections of 
relative order v^, we have 

{Of''{'S,)) = (2J+l)(Of«(^5o)>. (19) 

Identifying {O^"" Si)) = (0£), one finds that the corresponding term 
in the factorization formula ( |l6| ) reproduces the color-octet term in the 
nonrelativistic factorization formula (^5|) for P-wave production. 



The factorization formula (16) is not particularly useful in its most 
general form, since it involves infinitely many nonperturbative factors 
(Of). However, one can deduce from NRQCD how the various matrix 
elements scale with v for any particular quarkonium state H . For exam- 
ple, the scaling of {O^ (^^^^Lj)) with v is determined by the number 
of electric dipole and magnetic dipole transitions that are required to 
go from the dominant Fock state of the meson 77 to a state of the form 
|cc(n, ^'^+^Lj) -fgluons). The matrix element scales as ?;^+^^, multiplied 
by for each electric dipole transition and for each magnetic dipole 
transition. 

The relative importance of the various terms in the factorization 
formula ( |l6| ) is determined by the order in v of the matrix element (Of) 
and by the order in as(mc) of the parton cross sections da{cc{n) -\- X). 
The NRQCD factorization framework acquires predictive power when 
the expansion (|l^) is truncated to some low order in . If we truncate 
( p^ to lowest nontrivial order in v, regardless of the order in as, we 
recover the factorization formulas (^ for S-waves and ( p^ ) for P-waves. 
However, if the cross sections da for these terms are suppressed, other 
terms in the factorization formula may be important. This observation 
is the basis for a suggestion by Braaten and Fleming [M that color-octet 
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mechanisms may also play an important role in S-wave production. The 
leading color-octet matrix elements for the ■0 are (Og (^S'o)), (Og (^S'l)), 
and (Og (^Pj)), all of which are suppressed by relative to the leading 
color-singlet matrix element {Of {"^ Si)) . 

If the parton cross sections da multiplying the color-octet matrix ele- 
ments have the same magnitudes as the parton cross section multiplying 
(of then the corrections from the color-octet terms should be 

small. However, if the parton cross section multiplying the color-singlet 
matrix element is suppressed by powers of as{mc) or by small kinemati- 
cal parameters such as rn^/p^, the color-octet terms could dominate. An 
example is the gluon fragmentation function for producing a "0, which 
is discussed in Section 4.3. 

The implications of color-octet mechanisms for the production of 
quarkonium in high energy colliders are discussed in Section 5. There 
have also been many recent investigations of color-octet production 
mechanisms in lower energy experiments, including fixed target ttN and 



pN collisions |31[, _B-meson decays [52 , e+e annihilation |B3|, and pho- 




toproduction |54, Ka. These applications he outside the scope of this 



4.3 Fragmentation Functions for Heavy Quarkonium 

The NRQCD factorization formula implies that fragmentation func- 
tions for charmonium have the general form 

D,^H{z,^i) = ^d,^„(z,M) (Of), (20) 

n 

where di^n{z, /i) gives the probability for the parton i to form a jet that 
includes a cc pair in the state labeled by n, and (Of) is proportional 
to the probability for a pointlike cc pair in the state n to bind to form 
a charmonium state H . The coefficient di^„(z, 2mc) at the initial scale 
/i = 2mc involves only momenta of order TOcj and can therefore can be 
calculated as a perturbation expansion in the running coupling constant 
as{2mc). The fragmentation function can be evolved up to higher scales 
by using the evolution equations (0). Note that the only dependence 
on the quarkonium state H on the right side of ( po|) resides in the matrix 
elements. 

If the dominant Fock state for the meson H is |cc(l, ^'^^^Lj)), the 
fragmentation function in the color-singlet model is obtained by keep- 
ing only the term in (|2^) that involves (Of (^■^+^Lj)). In the case of 
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gluon fragmentation into Xcj, the coefBcient of {O^"'' i^Pj)) is logarith- 
mically infrared divergent at leading order in a^. The corresponding 
parton process is 5 — > cc(l_, ^Pj) + 5, and the divergence arises when 
the final state gluon becomes soft. In this region of phase space, the 
short-distance part of the fragmentation process is g ^ cc.{8,^ Si), and 
the divergence from the radiation of the soft gluon should be absorbed 
into the matrix element {0^'''{^Si)). The resulting expression for the 
fragmentation function at leading order in as is | |3C| ] 

+ ^^^^i^si^-^)(or'Cs,)), (21) 

where dj{z) is a dimcnsionlcss function of z. Since the matrix elements 
are the same order in v and the color-octet term is lower order in as , it 
can be expected to dominate. 

The process g cc{%,^Si) gives a term of order as in the gluon 
fragmentation function for any quarkonium state H: 



D,^H{z.2m^) « {O^i'Si)). (22) 

The delta function in (^2|) should be interpreted as a distribution in z 
that is peaked near z = \ with a width of order u^. All other terms in the 
fragmentation function have coefficients (ig_,„(2) of order or higher, 
or else have matrix elements that are higher order in than {O^ Si)) . 
The importance of the term (|2^) in the fragmentation function depends 
on the charmonium state H . Surprisingly, it may be the most important 
term in the fragmentation functions for the S-wave states ■(/' and if)' . 
Although the matrix element {O'^i^Si)) is suppressed by v"^ relative 
to {0^{^ Si)) , the leading term in the coefficient of the color-singlet 
matrix element comes from the parton process g — s- cc(X,^Si) + gg and 
is therefore of order a^ . The suppression of the color-singlet term by 
may be more effective than the suppression of the color-octet term by 
w**. The relative importance of the two terms can only be assessed after 
determining the value of the color-octet matrix element. An estimate 
for this matrix element will be obtained from Tevatron data in Section 
5.1. 
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5 COLLIDER APPLICATIONS 

Our understanding of heavy-quarkonium production in high energy col- 
Uders has been completely revolutionized by the theoretical develop- 
ments described in Sections 3 and 4. It is clear that fragmentation 
contributions can dominate at large transverse momentum, even if they 
are higher order in Ug- In addition, the formation of quarkonium from 
color-octet QQ pairs can in some cases dominate over the color-singlet 
contributions. In this Section, we discuss several applications to collider 
physics in which these developments play an important role. These ap- 
plications are the production of prompt charmonium at the Tevatron, 
the production of bottomonium and prompt charmonium at LEP, and 
the production of the Be in high energy colliders. 

5.1 Prompt Charmonium at Large px in pp Collisions 

Much of the impetus for the recent theoretical developments in heavy 
quarkonium production has come from Tevatron data on prompt char- 
monium production at large pr- As discussed in Section 2.3, the ex- 
perimentally measured rate at large px is orders of magnitude greater 
than the predictions of the color-singlet model at leading order in as- 
However, the Tevatron results can be explained by a combination of the 
fragmentation and color-octet mechanisms discussed in Sections 3 and 
4. 

In 1993, Braaten and Yuan |^ pointed out that fragmentation 
should be the most important charmonium production mechanism at 
sufficiently large pr- Doncheski, Fleming, and Mangano carried out 
the first explicit calculation of the fragmentation contribution to prompt 
charmonium production. They used the gluon and charm quark frag- 
mentation functions from the color-singlet model to calculate the frag- 
mentation contributions to direct ip and direct ip' production at large 
Pt at the Tevatron. They found that, although fragmentation does in- 
deed dominate over the leading-order color-singlet contribution for px 
greater than about 7 GeV, the predictions still fell more than an order 
of magnitude below the CDF data for inclusive ^ and ip' production. 

In subsequent calculations of prompt ip production | |5^ , contribu- 
tions from fragmentation into direct Xc followed by the radiative decay 
of Xcj into ip, were included. These contributions are roughly an order 
of magnitude larger than the direct ip contributions, and are dominated 
by the color-octet term in the gluon fragmentation function for XcJ- The 
value of the NRQCD matrix element {O^^^^^Si)) that was used in this 



QUARKONIUM PRODUCTION 



27 



calculation was estimated from data on B-mcson decays [||. The in- 
clusion of the contribution from gluon fragmentation into Xcj brings the 
theoretical prediction for inclusive ip production to within a factor of 3 
of the CDF data . Because of the many theoretical uncertainties that 
enter into the calculation, this factor of 3 discrepancy was considered 
acceptable at the time. However, in the case of prompt ■(/'' production, 
the theoretical prediction remained about a factor of 30 below the data. 
This dramatic discrepancy became known as the "CDF %p' anomaly" . 

The main difference between the prompt production of i]j and ij)' is 
that the ?/; signal is fed by direct Xcj's, while the ip' signal is not fed 
by any known charmonium states. The ■0' anomaly could be explained 
if there exist undiscovered charmonium states with sufficiently large 
branching fractions into ■0'. Among the possibilities that have been 
considered are D-wave states, higher P-wave states, and "hybrid" states 
whose dominant Fock state is ccg ^ . The main difficulty with these 
proposals is explaining why these states should have the large branching 
fractions into ip' that would be required to explain the data. 

An alternative solution to the -0' anomaly, based on the NRQCD 
factorization formalism, was proposed by Braaten and Fleming [50[ |. 
They suggested that the dominant contribution to the production of 
0' with large px comes from the color-octet term (^2|) in the gluon 
fragmentation function for -0'. This term represents the fragmentation 
of a gluon into a cc pair in a color-octet ^Si state at short distances, 
followed by the formation of the -0' from the cc pair through nonpertur- 
bative QCD interactions. The probability for the formation of the -0' is 
proportional to the NRQCD matrix element (Og (^^i)). While the nor- 
malization of this term in the cross section for tp' production depends 
on the undetermined matrix element (Og i^Si)), the pr-dependence 
of this contribution is predicted and found to be in good agreement 
with the CDF data By fitting to CDF data, one finds that the 

value of the matrix element is {of {^Si)) w 0.0042 GeV^ This value 
is consistent with the NRQCD prediction that it should be suppressed 
by v'^ with respect to the corresponding color-singlct matrix element 
{o f {^Si)) « 0.573 GeV^ Thus this solution to the 0' anomaly is at 
least plausible. 

This proposal was given further support by subsequent CDF data 
in which the prompt '0's were separated into those that come from Xc 
decay and those that do not |l^]. The observed rate for prompt V''s 
not from Xc decay was about a factor of 30 larger than the theoretical 
predictions for direct ip from fragmentation in the color-singlet model. 
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Cacciari, Greco, Mangano, and Petrelli showed that this data could be 
explained by color-octet gluon fragmentation As in the case 

of ip' , the p^-dependence predicted by color-octet fragmentation is in 
good agreement with the data. The value of the color-octet matrix 
element that is required to fit the data is (0^(^51)) w 0.014 GeV^. 
This is small compared to the corresponding color-singlet matrix element 
(of (35'i)) « 1.00 GeV^ consistent with the suppression by predicted 
by NRQCD. 

The simple picture that seems to emerge from the CDF data is that 
the dominant mechanism for the production of "0, ^' , or Xcj at large pt 
is the pcrturbative fragmentation of a gluon into a cc pair in a color-octet 
'^Si state, followed by the nonperturbative formation of charmonium 
from the cc pair. The differential cross section therefore reduces to that 
for producing a gluon convoluted with a fragmentation function: 

da{pp ^ H{pt) + X) = y2 dxifj/p{xi) dx2 fk/p{x2) 

Jk -^0 -^0 

X dz dd{jk g{pT/z) + X, ^ft-ag) Dg^niz, Aiftag)- (23) 
Jo 

The dominant term in the gluon fragmentation function is assumed to 
be the term proportional to {Og {'^Si)). It is given by (^2|) at the scale 
/^frag = 2mc, but it should be evolved up to the scale /ifrag — Pt/z 
using the evolution equations (0). Of course, a thorough analysis of 
charmonium production at large pt must include additional terms in 
the gluon fragmentation function that are of order aj. and higher. It 
must also include the contributions from the fragmentation of other 
partons, such as the charm quark and light quarks. 

In Figures 1,2, and 3, the predictions of this simple picture of char- 
monium production are compared with the most recent CDF data on 
prompt ip not from Xc decay ||l^ , prompt ?/' from Xci and Xc2 decay |l^ , 
and prompt -0' respectively. The predictions of the color-singlet 
model at leading order in ag are shown as dashed lines and they fall 
orders of magnitude below the data. The dotted lines in Figures 1 
and 2 are the fragmentation contributions calculated in the color-singlet 
model. While they increase the theoretical predictions by more than an 
order of magnitude at the largest values oipx: they still fall about a 
factor of 30 below the experimental measurements. The solid curves in 
Figures 1,2, and 3 are the contributions from color-octet fragmentation, 
with the color-octet matrix elements {O^ {^Si)) adjusted to fit the data. 
The solid curves differ only in their normalizations. The normalizations 



QUARKONIUM PRODUCTION 



29 



in Figures 1, 2, and 3 are proportional to {Ofi^Si)) + 0.57(of (^S'l)), 
0.27(08^- +0.14(O|-(35i)) = 1.51(08^-(35i)), and (ofi'S^)), 
respectively. Their shapes agree reasonably well with the data. The 
values of the matrix elements for and ip' are given above. The value 
{O^'^i^Si)) = 0.0076 GeV^ obtained for the Xc matrix element is 
significantly larger than the most recent value obtained from i?-meson 
decay ijEoj . In calculating these curves, we imposed a pseudorapidity cut 
of I?;! < 0.6 on the charmonium state and we used the MRSDO parton 
distribution set. The scales /i^j, fip, and iJ-frag were all chosen to be 
equal to the transverse momentum pr/z of the fragmenting giuon. 

While the fragmentation contribution must dominate at sufficiently 
large pr, the contributions to the cross section that are suppressed by 
factors of ra^/p^ may be important at the values of pt that are experi- 
mentally accessible. These contributions have been included in a recent 
calculation by Cho and Leibovich Q. They calculated the cross sec- 
tions for all 2 — > 3 parton processes of the form ij —^cc+k that produce 
cc pairs in color-octet states with angular momentum quantum numbers 
•^5*1, ^5*0, and '^Pj. They used these to calculate the cross section for di- 
rect ip and ip' production, including all color-octet contributions that are 
suppressed only by w^. By fitting the pr distributions for all pt greater 
than 5 GeV, they extracted values for the matrix elements {Of(}So)) 
and {Ofi^Pj)), as well as {0^{^Si)). The best fit requires a value 
for that is significantly smaller than that obtained above. 

The terms involving {Og(^So)) and {Of^^Pj)) dominate in the range 
5 GeV < PT < 10 GeV. At pr > 10 GeV, the term involving (Oti^Si)) 
dominates and it asymptotically approaches the contribution from the 
leading color-octet term in the gluon fragmentation function. The dif- 
ference between the full calculation of the differential cross section and 
the fragmentation approximation is less than 20% at pT = 10 GeV. 

The greatest weakness of the proposal to explain the CDF data by 
color-octet production mechanisms is that the color-octet matrix ele- 
ments are free parameters that can be adjusted to fit the data. The 
predictive power in this proposal lies in the fact that the NRQCD ma- 
trix elements are universal, and will also enter into other charmonium 
production processes. The most convincing proof that the color-octet 
proposal is correct would be to show that the values of the matrix ele- 
ments determined from the CDF data are required to explain charmo- 
nium production in other high energy processes. One possibility is 
decay, which will be discussed in Section 5.2. 

The color-octet fragmentation mechanism also gives other predic- 
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tions that can be tested experimentally. The most dramatic prediction 
is that direct ip^s and ip'^s at large px will have a large spin alignment. 
Cho and Wise ||6^ pointed out that ip's and ip'^s produced by color- 
octet gluon fragmentation tend to inherit the transverse polarization of 
the nearly on-shell gluon. The process g cc{8,^Si) produces a cc 
pair whose total spin is transversely polarized. The approximate heavy- 
quark spin symmetry of NRQCD implies that the spin state of the cc 
pair will be affected very little by the nonperturbative QCD effects that 
are involved in binding the cc pair into a ■0. At leading order in as, the 
predicted polarization is 100% transverse. Beneke and Rothstein an- 
alyzed the radiative corrections and concluded that they can decrease 
the polarization by only about 10% j6^. Cho and Leibovich have 
shown that the dominant corrections to the spin alignment at Tevatron 
energies come from terms in the cross section that involve (Og (^^o)) 
and (Oti^Pj)) and fall as ni^/p^. These corrections remain small at 
the largest values of px measured at the Tevatron. A large transverse 
spin alignment is therefore a robust signature for the dominance by 
color-octet fragmentation of the production of direct i/j's and ip"s at 
high pt- 

There are other predictions of the color-octet fragmentation mech- 
anism that can be tested at the Tevatron. As pointed out by Barger, 
Fleming, and Phillips, the production of two quarkonium states with 
large pt will be dominated by the production of two high-p^ giuons 
which both fragment into quarkonium [|4|. Taking into account the 
color-octet mechanism, they found that ipip events at large px should 
be detectable at the Tevatron. The color-octet fragmentation mecha- 
nism also gives predictions for correlations between charmonium at large 
Pt and the other jets produced by the pp collision, but these have not 
yet been explored. 

The differential cross sections for the production of bottomonium 
states at the Tevatron have also been calculated by Cho and Leibovich 
pO| , |6]| . The predictions of the color-singlet model fall more than an 
order of magnitude below the CDF data for px > 5 GeV . The pro- 
duction rates can, however, be explained by including color-octet terms 
in the cross section for the states T{nS) and Xbj{n), n — 1,2,3. The 
values of the color-octet matrix elements that are required to fit the 
CDF data are consistent with expectations from NRQCD. The calcula- 
tions of Cho and Leibovich show that fragmentation is not important 
in bottomonium production at the Tevatron. The error from neglecting 
contributions that fall as m^/p^ is greater than 100% even if pT is as 
large as 20 GeV. 
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Color-singlet and color-octet fragmentation processes in inelastic 
ip photoproduction at the HERA ep collider have also been investi- 
gated The dominant contributions for values of pr accessible 
to experimental studies appear to be color-singlet charm quark frag- 
mentation and color-octet processes that are suppressed by m^/p^. In 
the color-octet contributions the cc pair is created at short distances in 
either a ^5*0 or ^Pj state. Because the numerical values of the color- 
octet matrix elements (Og (^S'o)) and {Of{^Pj)) are not measured very 
precisely it is not possible to determine which process is the dominant 
one at HERA. However at asymptotic values of pr the fragmentation 
contribution will be the most important contribution. Thus for large 
Pt data from HERA can be used to study a fragmentation mechanism 
that is not so easily probed at the Tevatron and at LHC. 

5.2 Charmonium and Bottomonium in Decay 

The decay of the provides a laboratory for the study of heavy quarko- 
nium with large transverse momentum that is complementary to pp col- 
lisions. The color-singlet model predicts that the production rates for 
prompt charmonium and for bottomonium are so small that they are 
unlikely to be observed at LEP. However, the dramatic failure of the 
color-singlet model in pp collisions suggests that its predictions for 
decay should also be reexamined. 

The predictions of the color-singiet model for charmonium produc- 
tion in decay have been studied thoroughly |||, |7[ |6^ . The largest 
production process for in the color-singlet model is Z^ — > ipcc. The 
rate for this process was first calculated in 1990 by Barger, Cheung, 
and Keung p^ , who found it to be surprisingly large. It is two orders 
of magnitude larger than the rate for Z*^ ipgg [Q, despite the fact 
that both processes are of the same order in as- This was explained in 
1993 by Braaten, Cheung, and Yuan, who pointed out that the reason 
Z'^ —> Ipcc is so large is that this process has a fragmentation contri- 
bution that is not suppressed by m^/M| They showed that the 
dominant terms can be factored into the decay rate for Z'^ — > cc and 
fragmentation functions for the processes c ipc and c — > tpc'- 

^(ZO ^ +X) ^ 2 T{Z" ^ cc) D,^^{z^), (24) 

where — 2E^/Mz- Since the momentum scale for the fragmenta- 
tion process is of order mc, the rate is proportional to alirric) instead 
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of a^{Mz)- The resulting prediction for the branching fraction for di- 
rect tfj production in Z'^ decay is 2.9 x 10^^. The experimental study of 
prompt charmonium in Z° decay is complicated by the large background 
from the decay of 5-hadrons, which must be removed using vertex detec- 
tors. There are preliminary results from LEP | ]69| that indicate that the 
branching fraction for prompt ip production is around 10^^, well above 
the predictions of the color-singlet model. 

Another fragmentation process in the color-singlet model was con- 
sidered by Hagiwara, Martin, and Stirling in 1991 They calculated 
the rates for the process —> qqg*, with charmonium being produced 
by the virtual gluon through the processes g* tp99 and g* Xcjg- 
The rate is at least an order of magnitude smaller than the charm frag- 
mentation process described above, but the experimental backgrounds 
are less severe. The DELPHI collaboration has recently reported 
4.1 X 10^"* as a limit on the branching fraction for — > ipX from this 
process. 

The color-octet mechanism that was introduced to explain prompt 
charmonium production at the Tevatron also offers new possibilities for 
charmonium production at the Z'^ resonance. Cheung, Keung, and Yuan 
and Cho [j7l| have studied prompt charmonium production at LEP 
via the color-octet mechanism. The color-octet process that is leading 
order in as involves the short-distance decay Z° — > cc(8, '^Si)+g, but it 
has a negligible branching fraction because it is suppressed by a short- 
distance factor of m^/M|. The dominant color-octet process is of order 
and involves the short-distance decay Z^ cc{8,^Si) + qq, where q 
can be one of the light quarks u, d, or s or one of the heavy quarks c 
or b. Apart from different coupling constants, color factors, and long- 
distance matrix elements, this process is identical to the electromagnetic 
production of ip studied by Fleming and discussed in Section 3.2. 
At leading order in as, the perturbative QCD factorization formula for 
this process has the form 

—{z°^ij{z^) + x) « 2^r(z°^gg)D,^v(^^:M) 

azTh ^ 

d2/^(^° ^ .9(z^/j/)gg,Ai)i^<,^^(2/), (25) 

where — 2E^/Mz and Zg = ZTp/y. The fragmentation function 
Dg^jl,{y) for the formation of a V' in a gluon jet through the color-octet 
mechanism is given in (p^), except that it must be evolved up to the scale 
E^/y set by the gluon energy. The fragmentation function Dq^^{z) for 
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the formation of a V' in a light quark jet can be obtained from (|l3|) by 
replacing 0^3% by aHpt{^ S^)) I 

In Figure 4, we compare the ■0 energy distribution predicted by the 
color-singlet model with the distribution from the most important color- 
octet process Z'^ cc(8, ^Si)+qq. The most striking feature in Figure 4 
is that the energy distribution of the color-octet process dominates over 
that of the color-singlet model for all values of = 2i?^ /AIz ■ There is 
also a dramatic difference between the shapes of the two distributions. 
The energy distribution from the color-singlet process is rather hard 
because of the nature of heavy quark fragmentation. The distribution 
for the color-octet process is very soft and has a pronounced peak near 

= 0.1. The prediction for the branching ratio ior Z'^ ^ tp + X from 
the color-octet process is 1.4 x 10~^. This is consistent with the DELPHI 
bound for the color-singlet process Z^ — > tpggqq, which produces V''s 
with a similar event topology. 

The color-octet fragmentation process is also important in the pro- 
duction of V'' and Xc- The energy distributions are predicted to be 
similar to that for the "0 shown in figure 4. The production rates are 
predicted to be smaller by about a factor of 3 for the ip' and larger by 
about a factor of 5 for the three spin states of the Xcj combined. 

Just as in the case of pp annihilation, the color-octet fragmentation 
mechanism makes definite predictions for the spin alignment of the tp 
in Z° decay. The color-octet process Z° cc{8,^Si) + qq involves 
contributions from both gluon fragmentation and light quark fragmen- 
tation. At leading order in as, gluon fragmentation only produces tp^s 
that are transversely polarized. However, a light quark can fragment 
into a longitudinally polarized ^ even at leading order in as . This leads 
to significant degradation of the spin alignment of the ip. The spin align- 
ment varies significantly with z, but the average transverse polarization 
may be as small as 75% fz^ . 

Bottomonium states can also be produced in Z° decay. These states 
are simpler to study experimentally, since there is no background to bot- 
tomonium production from non-prompt production mechanisms analo- 
gous to 6-hadrons decay into charmonium. A signal from the J^*^ — 1 
states T(15'), T(2S'), and T(35') has in fact been observed at LEP. The 
preliminary result from the OPAL detector is that the branching frac- 
tion summed over the three states is | |7^ 

3 

Br(Z" ^ T{nS) + X) = il.2toi ± 0.2) x 10"^ . (26) 

n=l 
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The direct production mechanisms for T(nS) in decay are essentially 
identical to those for prompt -0. The dominant color-singlet process is 
the decay Z° — > T{nS) + bb. This process has a fragmentation contribu- 
tion which corresponds to the decay Z° bb followed by the formation 
of T{nS) through the fragmentation of the b or 6. The prediction for 
the sum of branching fraction given in ( |2^ ) is 1.6 x 10~^, which is about 
an order of magnitude below the OPAL result. Color-octet mechanisms 
for the production of T{nS) and Xbjin) have been studied by Cho Q. 
The dominant mechanism is the formation of bottomonium through 
the short-distance process Z° bb{8,^Si) + qq, q — u,d,s, which in- 
volves both gluon fragmentation and light quark fragmentation. Using 
values for the color-octet matrix elements (O^j^Si)) for the bottomo- 
nium states that were extracted from CDF data [£5|, Cho predicted 
that the branching fraction in ( p6|) should be 4.1 x 10^^. Given the 
large uncertainties in the color-octet matrix elements, this is in reason- 
able agreement with the OPAL result (p6|). If one can show that the 
same values of the color-octet matrix elements explain T production in 
pp colliders and in Z° decay, this would provide strong support for the 
proposal that the production is dominated by color-octet mechanisms. 

5.3 Production of the Be 

In the Standard Model, the only bound states consisting of two heavy 
quarks with different flavors are the be mesons. The ground state of the 
be system is the pseudoscalar state Be- The mass spectrum of the be 
mesons can be predicted reliably from quark potential models that are 
tuned to reproduce the spectra of charmonium and bottomonium Q . 
According to the potential- model calculations, the first two sets of S- 
wave states, the first and probably the entire second set of P-wave states, 
and the first set of D-wave states all lie below the BD flavor threshold. 
Since QCD interactions are diagonal in flavors, the annihilation of be 
mesons can only occur through a virtual and is therefore suppressed 
relative to the electromagnetic and hadronic transitions to lower-lying be 
states. Thus all the excited states below the BD threshold will cascade 
down to the ground state Be via emission of photons and/or pions. A 
calculation of the inclusive production of the B^ meson must therefore 
include the contributions from production of all the S-wave, P-wave, 
and D-wave states that are below the BD threshold. 

The production of the Be and the lowest '^Si state B* in e^e~ an- 
nihilation was first computed to leading order in as in the color-singlet 
model by Clavelli and by Amiri and Ji iTq] . The lowest-order process for 
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Be production is e+e^ Bc + bc. Chang and Chen |^ showed that the 
energy distributions for direct B^. and B* production could be expressed 
in terms of the fragmentation functions Di^g^{z) and Di^gt {z). For 
example, the production rate for Be at the resonance in the scaling 
limit Mz oo with z = 2Eb^/Mz fixed has the form 

^{Z^ ^ Be{z)+X) « r(Z° ^ 66) D-,^B^[z). (27) 

The fragmentation functions D^^g (z) and Di_^g,(z) were also calcu- 
lated later by Braaten, Cheung, and Yuan who pointed out that 
they were universal and could also be used to calculate direct Be pro- 
duction in other high-energy processes, such as pp collisions. The first 
calculations of Be production in hadron colliders using fragmentation 
functions were carried out by Cheung j?^ . The fragmentation func- 
tions were subsequently calculated for the P-wave states of the be system 
p^ , |3^ and even for the D- waves |36|, allowing a complete calculation 
of the inclusive Be production rate. The advantage of the fragmentation 
approach over carrying out a full calculation 

of the production rate to leading order in is that the calculation 
of the fragmentation function is much simpler and it can be extended 
to P-wave and D-wave states relatively easily. The limitation of the 
fragmentation approach is that it is accurate only at sufRciently large 
Pt- 

Color-octet production mechanisms are expected to be much less 
important for be mesons than for charmonium or bottomonium. All 
terms in the gluon fragmentation functions are of order or higher. 
For the 6-quark fragmentation functions, the leading color-singlet terms 
and the leading color-octet terms are both of order a^. For S-wave 
states like Be and B*, the color-octet terms are therefore suppressed 
relative to the color-singlet terms by a factor of from the ratio of 
the NRQCD matrix elements. For P-wave states, a color-octet term in 
the fragmentation function is necessary to avoid an infrared divergence 
in the color-singiet fragmentation function. However, this divergence 
occurs first at ncxt-to-leading order in as- Thus the color-octet term 
is not as important as it is for the case of gluon fragmentation into Xc, 
where the divergence occurs at leading order. 

The production of the Be and B* at hadronic colliders like the Teva- 
tron and the Large Hadron Collider (LHC) has recently been calculated 
by several groups, both through complete C(af) calculations |7^ 
and using the simpler fragmentation approximation ||7^, ^ . Discrepan- 
cies among the earlier calculations raised questions about the accuracy 
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of the fragmentation approximation in the be case. Recent calculations 
have shown that the fragmentation approximation is very accurate 
for the Be for pj- as low as 10 GeV. However, for the B*, the accuracy 
is much lower, decreasing to 20% only for around 30 GeV. For 
production in 77 collisions, even larger values of pr are required be- 
fore the fragmentation approximation becomes accurate |Q. The fact 
that such large values of px are required in order for fragmentation 
to dominate can be attributed to the presence of an additional small 
parameter nic/mb in the problem. Contributions that fall off asymptot- 
ically as /p^ may dominate at sub-asymptotic values of pt if they 
are enhanced by powers of mi,/mc. 

The only calculations presently available for the production of the 
P-wave be mesons have been carried out using the fragmentation ap- 
proximation. This approach has been used by Cheung and Yuan to 
calculate the inclusive production of the Be in pp collisions, including 
the contributions from all the S-wave and P-wave states below the BD 
threshold |7^. With acceptance cuts of pT > 6 GeV and \y\ < 1, the 
inclusive production cross section for the B^ meson at the Tevatron is 
about 5 nb. The contributions from the excited S-wave states and from 
the P-wave states are 58% and 23%, respectively, and hence significant. 
The D-wave states were not included in this analysis since they are ex- 
pected to contribute only about 2% to the inclusive production of B^ 
. The search for the Be is now underway at the Tevatron, and pre- 
liminary results have already been presented by the CDF collaboration 
If the Be is not found in the present run of the Tevatron, it should 
certainly be discovered after the installation of the Main Injector, which 
will boost the luminosity by about a factor of 10. 

The Be may also be discovered at LEP. Candidate events for the 
Be have already been reported by the ALEPH collaboration |Q. The 
branching fraction predicted from the color-singlet processes Z'^ ^ Bc + 
be and Z° ^ B* + be are 5.9 x lO"'^ and 8.3 x 10"^, respectively 
If the fragmentation functions are used to estimate the contributions 
of higher S-wave, P-wave, and D-wave states, the inclusive branching 
fraction for Be is predicted to be larger by about a factor of 5. 

6 OUTLOOK 

The NRQCD factorization formalism provides a very general framework 
for analyzing the production of heavy quarkonium. It implies that the 
production process is much more complex than has been assumed in the 
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color-singlet model. There are infinitely many nonperturbative matrix 
elements that contribute to the cross sections, although only a finite 
number of them contribute at any given order in . The most dramatic 
consequence of this formalism is that color-octet mechanisms sometimes 
give the largest contributions to the cross section. In the case of P- waves, 
there is a color-octet term at leading order in that must be included 
for perturbative consistency. In the case of S- waves, color-octet terms 
are suppressed by w^, but they may be important if the color-singlet 
term is suppressed by other small parameters such as a., . 

The NRQCD factorization formalism suggests an explanation for 
the large production rates for bottomonium and prompt charmonium 
that have been observed at the Tevatron. They can be attributed to 
color-octet terms in the cross section that are important because the 
color-singlet terms are suppressed by powers of and mQ/p^. The 
magnitudes of the color-octet matrix elements that are required to ex- 
plain the Tevatron data are in accord with expectations from NRQCD. 
This explanation leads to many predictions that can be tested exper- 
imentally. The ultimate test is that the same matrix elements must 
be able to explain heavy quarkonium production in other high energy 
processes, such as Z'^ decay. 

In the NRQCD factorization formalism, there are many nonpertur- 
bative matrix elements that must be determined phenomenologically in 
order to make theoretical predictions for quarkonium production. For- 
tunately, there is a wealth of experimental data that can be used to de- 
termine these matrix elements. Quarkonium is produced as a byproduct 
of almost every high energy experiment. In this review, we have only 
discussed data from the highest-energy colliders, the Tevatron and LEP. 
However there is also data on quarkonium production from lower-energy 
e+e~ and ep colliders. There is also an abundance of data on pN, ttN, 
and "/TV collisions from fixed target experiments. By carrying out a com- 
prehensive analysis of all the data on quarkonium production, it should 
be possible to determine most of the NRQCD matrix elements that are 
of phenomenological importance. 

An enormous amount of theoretical work will be required in order to 
carry this program to completion. For every process, the production rate 
should be calculated to next-to-leading order in all the small parameters 
in the problem, including v^, as, and kinematical parameters such as 
^q/Pt- Most calculations until very recently were carried out within 
the color-singlet model, and therefore include only contributions that 
are of leading order in v^. Even in the color-singlet model, the only 
process for which the corrections at next-to-leading order in have 
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been calculated is photoproduction of J/'0 The effort to calculate 
contributions that are beyond leading order in u^, such as color-octet 
production mechanisms, has just begun. 

While the NRQCD factorization approach provides a very general 
framework for analyzing quarkonium production, it should not be re- 
garded as a complete theory. The derivation of the factorization for- 
mula breaks down when there are hadrons in the initial state if the 
quarkonium is produced with small transverse momentum. Therefore, 
processes involving diffractive scattering, such as the elastic scattering 
process 7p — > ipp, are not described by this formalism. The NRQCD 
formalism also is not appropriate for describing the formation of char- 
monium from intrinsic cc pairs that come from the parton distribution 
of a colliding hadron iQ . 

The complications of diffractive scattering and intrinsic cc pairs do 
not arise in the production of heavy quarkonium at large px , and there- 
fore the theoretical analysis of this process is particularly clean. The 
factorization theorems of perturbative QCD guarantee that the produc- 
tion is dominated by fragmentation, the formation of heavy quarkonium 
within "jets" that are initiated by single high-energy partons. The large 
cross sections for prompt charmonium production that have been ob- 
served at the Tevatron can be explained by including a color-octet term 
in the gluon fragmentation function. While the normalization of the 
cross section is not predicted, this mechanism does give other predic- 
tions that can be tested experimentally. The most dramatic prediction 
is a large transverse spin alignment for direct -0 and '0' at large px, 
but the spin alignment has not yet been measured. This mechanism 
also gives predictions for the correlations between high-p^ charmonium 
states and other jets produced by the collision that can be tested. 

An enormous amount of theoretical work remains to be done in or- 
der to obtain precise predictions for fragmentation that can be compared 
with experiment. Quantitative estimates of all the relevant NRQCD ma- 
trix elements are required in order to determine which terms in the frag- 
mentation functions are numerically important. The important terms 
should all be calculated to next-to- leading order in as . It is also impor- 
tant to calculate the power corrections to the cross sections that fall as 
ruQ/p^, since they may give the largest corrections to some observables. 

New experimental data from high-energy colliders is providing strin- 
gent tests of our understanding of the production of heavy quarkonium. 
The NRQCD factorization approach suggests that cross sections should 
be calculable in terms of the heavy-quark mass mg, the running cou- 
pling constant and a few matrix elements that can be determined 
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phcnomenologically. This approach provides an explanation for the large 
production rates that have been observed at the Tevatron and at LEP, 
but further effort, both theoretical and experimental, will be required 
to show conclusively that this explanation is correct. If this effort is 
successful, it will mark a major milestone on the road to a compre- 
hensive description of heavy-quarkonium production in all high energy 
processes. 
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Figure Captions 

Figure 1 CDF data on the differential cross section for prompt ^'s 
that do not come from Xc decay as a function of pr- The curves are 
the leading-order predictions of the color-singlet model (dashed curve), 
the predictions from fragmentation in the color-singlet model (dotted 
curve), and the contribution from gluon fragmentation via the color- 
octet mechanism (solid curve) with the normalization adjusted to fit 
the CDF data. 

Figure 2 CDF data on the differential cross section for prompt ip^s 
from Xc decay as a function ofpx- The dashed and solid curves are as 
described in Figure 1. 

Figure 3 CDF data on the differential cross section for prompt ip'^s as 
a function of pt- The curves are as described in Figure 1. 

Figure 4 Predictions for the energy spectrum dV/dz for prompt ■(/''s 
from decay in the color-singlet model (dashed line) and from the 
leading color-octet process (solid line). 
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